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A nonlinear beatwave regime of plasma wave excitation is considered. Two beatwave drivers are 
considered: intensity-modulated laser pulse and density-modulated (microbunched) electron beam. 
It is shown that a long beatwave pulse can excite strong plasma waves in its wake even when 
the beatwave frequency is detuned from the electron plasma frequency. The wake is caused by the 
dynamic bi-stability of the nonlinear plasma wave if the beatwave amplitude exceeds the analytically 
calculated threshold. In the context of a microbunched beam driven plasma wakefield accelerator, 
this excitation regime can be applied to developing a femtosecond electron injector. 

PACS numbers: 52.35.Mw, 52.38.Kd, 41.75.Jv 

Beatwave excitation of electron plasma waves continues attracting significant attention as a basic nonlinear plasma 
phenomenon, and as a viable approach to plasma-based particle acceleration 1, 2, 3, 4]. Beatwave excitation mech- 
anism is realized when the driver intensity (laser or particle beam) is modulated with the temporal periodicity 
of the plasma wave. The linear one-dimensional theory of the beatwave-driven plasma wave generation is well 
understood Q, 0|, and its most important predictions are as follows. First, the effectiveness of plasma wave ex- 
citation is strongly dependent on the difference between the beatwave frequency ujb and plasma wave frequency 
Wp = a/ A-Ke^riQ/m (where — e and m are the electron charge and mass, and no is the plasma density): the smaller is 
the frequency detuning Atj = lob — ujp, the larger is the resulting plasma wave inside the beatwave. Second, only if 
the beatwave pulse is short enough for its bandwidth to be comparable to Aw, an appreciable plasma wave is left in 
its wake. 

In this Letter I demonstrate that these conclusions are no longer valid when the relativistic nonlinearity of a plasma 
wave is accounted for. In particular, a strong plasma wave can be excited in the wake of a relatively long beatwave 
pulse of duration ti, S> 1/ Aw due to the nonlinear phenomenon of dynamic relativistic bi-stability (RB) ^] . Another 
manifestation of RB is that, at a certain critical strength of the beatwave driver, a weak driven plasma wave becomes 
unstable, and a much higher amplitude wave is excited. Linear estimates of the plasma wave amplitude fail when 
the beatwave amplitude exceeds this detuning-dependent critical strength. As the time-dependent beatwave strength 
increases and exceeds the critical value, significant pulsations of the plasma wave amplitude occur. These pulsations 
indicate that significant energy exchange takes place between the plasma wave and the driver. This effect can be 
exploited in a plasma wakefield accelerator driven by a microbunched electron beam 0]: bunches in the head of the 
beam excite while those in the back deplete plasma waves, thereby gaining energy. 

Relativistic bi-stability was originally described (3] for a magnetized electron subjected to cyclotron heating. Appli- 
cations of RB to electron cyclotron heating of fusion plasmas la Bl ha ve been later suggested. Although the nonlinear 
nature of electron plasma waves has been noted before [TllTlTiarr^. the RB of plasma waves has not been explored, 
either as a basic phenomenon or in the context of plasma-based accelerators. 

The one-dimensional relativistic dynamics of the cold plasma driven by a beatwave can be described using a 
Lagrangian displacement of the plasma element originally located at Zg: z{t) = Zq + (^{t,zo). It is assumed that 
the beatwave generated by either a pair of frequency-detuned laser beams, or a modulated electron beam, is moving 
with the speed close to the speed of light c, and, therefore, all beatwave quantities are functions of the co-moving 
coordinate r' = ujp{t — z/c) = t — lOpC^jc. Introducing the normalized displacement C ~ ^pC,lc and longitudinal 
relativistic momentum p = jd(/dT, where 7 = ^1 — iP/c^, equations of motion take on the form 

f = ^I=, f^-C + «(.') COS..'. (1) 

Assuming that |Aw| <^ ujp (near-resonance excitation), transverse momentum of the plasma has been neglected and 
the relativistic 7-factor simplified to 7 = ^1 + p^. The first term in the force equation is the restoring force of the 
ions, and the second term signifies the beatwave with the frequency lob = Loujp. The nonlinear in C modification of 
the beatwave in the rhs of Eqs. ^ is neglected in what follows. For a pair of linearly polarized laser pulses with 
electric field amplitudes Ei and E2 and the corresponding frequencies wi and 0^2 = loi — lob the normalized beatwave 
amplitude a = {e/mcYEiE2/2ujiUj2 ^3- ^'^'^ ^ driving electron bunch with the density profile rn, = nto + dnbSinujT 
it can be shown that a = Snb/no. Although arbitrary profiles of a(r) are allowed, it is assumed that \da/dT\ <^ \a\. 
The total energy density of the plasma wave Up/riomc^ = ^/l +p'^ + is changed via the interaction with the 
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FIG. 1: Steady-state solutions of a driven plasma wave as a function of the beatwave amplitude a. Solid lines 1, 2; stable 
equilibria for lo = 0.95; dashed line: unstable equilibrium for ui = 0.95; dot-dashed line: resonant excitation with ui — 1. 



beatwave. The effect of the plasma wave on the beatwave is neglected for the moment and addressed towards the end 
of the Letter. 

Although Eqs. Q can be solved numerically at this point, further simplification is made by assuming p — 
u cos {lut + 4>), where u and 4> are slowly varying functions of r. In the weakly relativistic approximation <C 1 
obtain: 

du a , . , 

- = -cos</> (2) 

u'^^-^sin^-^{uo^ -l + Zu'/S). (3) 
or 2 2uj 

Depending on the beatwave frequency oj and the amplitude a, equilibrium solutions du/ dr ~ (steady amplitude) 
and d(j)/dT = (phase-locking to the beatwave) of Eqs. (|2I3I) can have one or three real roots. For any uj there 
is a stable equilibrium point: (f)o — — 7r/2 and uq > found as the root of the third-order polynomial equation 
V(uo) = uo(w^ — 1 + 3/8mq) — Lua. For the most interesting a; < 1 regime additional solutions 0o = 7r/2 and uq > 0, 
where Uq is the positive root of V (uq) — —Lua, may be found, depending on the beatwave amplitude. Specifically, 
there are no additional positive roots for a > acrit, where Qcrit = 4-\/2(l — w^)^/^/9a;, and two positive roots ui,2 
for a < acrit (one of them unstable). Stable equilibrium amplitudes mq with 4>q = 7r/2 (Branch 1) and = — 7r/2 
(Branch 3), as well as the unstable one (Branch 2) are plotted in Fig. ^ as a function of the beatwave strength o for 
UJ = 0.95 {ami. = 0.02). Equilibrium bi-stability corresponding to Branches 1 and 3 is universal for any nonlinear 
pendulum |a. Il4|. including a weakly damped one. Equilibrium solutions are meaningful only if the plasma wave is 
phase- locked to the beatwave: d(p/dT « 0. As shown below, this is not the case when the peak beatwave amplitude 
exceeds acrit- Nonetheless, a dynamic RB described below occurs even in the absence of phase-locking. 

Consider plasma response to a Gaussian beatwave pulse a(T) — ao exp (— t^/t|), where tl » 1/|1 — cj] is the 
normalized pulse duration. For ao < acrit the plasma response is as follows: amplitude u adiabatically follows a(T) 
by staying on the Branch 1 and following the equilibrium trajectory schematically shown by arrows in Fig. ^ The 
adiabaticity condition is ^btl S> 1, where fls is the bounce frequency around the equilibrium point uq such that 
V{uq) = —ujair). Linearizing Eqs. H2I3|I around (j) = 7r/2 and u — uq yields fi^ — a{u1^^f. — Mq)/4cjmo, where 
Ucrit — 2^2(1 — 1^^)/3 is the critical plasma wave amplitude corresponding to the merging point between Branches 
1 and 2 in Fig. For oq < acrit plasma oscillation is indeed phase-locked to the beatwave at (/iq w 7r/2 during the 
ramp-up and most of the ramp-down of the laser pulse (although phase-locking is lost when the pulse amplitude 
becomes very small on the down-ramp). As the result, plasma wave amplitude returns to a very small value in the 
wake of the beatwave, as shown by a dot-dashed line in Fig. [3 The longer is the beatwave pulse duration t^, the 
smaller is the wake because its non- vanishing amplitude is due to the adiabaticity violation for finite . 

Situation changes for ao > acrit- as air) approaches acrit-, the adiabatic condition is violated (noted in the context 
of electron cyclotron heating |3), and phase-locking at <^o = ''■/2 is no longer possible. Thus, the transfer to 
Branch 3 schematically shown by a vertical arrow in Fig. ^ becomes feasible, and the plasma wave amplitude can 
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FIG. 2: Excitation of a plasma wave by a Gaussian heatwave pulse (dotted line), a(r) — ao exp [— r /tI], tl = 150. Solid 
line: ui = 0.95, above-threshold excitation with ao = 0.023 > acrit ~ 0.02; dashed line: resonant excitation with cu = 1 and 
ao — 0.023; dot-dashed line: u} — 0.95, below-threshold excitation with ao — 0.018 < acrit- 



dramatically increase. In the presence of a finite plasma wave damping this indeed happens: the subsequent decrease 
of the beatwave amplitude results in phase-locking at (f>o — — 7r/2, with u following along the Branch 3. Without 
damping, there is no mechanism for the plasma wave to reach the equilibrium amplitude given by the upper Branch 
3. As shown below, a conservation law prohibits the jump between Branches 1 and 3. 

Nevertheless, even without damping, a significant plasma wave is left behind the finite-duration beatwave pulse 
(Fig. [3 solid line). The previously unaccessible finite-amplitude solution has been reached due to the effect of the 
dynamic RB which is best understood through the conservation of the effective Hamiltonian of the driven plasma 
wave. The effective Hamiltonian 
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can be used to express Eqs. (I2I3|I in the form of u = {l/ujdH/dip, cj) = —{l/u)dH/du. For a slowly changing beatwave 
amplitude a(r) the Hamiltonian is almost conserved: dH/dr = Q.busaKpda/ dr w 0. This constitutes the conservation 
law preventing the jump between Branches 1 and 3. For the initially quiescent plasma a = and m = before the 
arrival of the beatwave. Therefore, H after its passage, as confirmed by numerical simulations of various pulse 
durations and amplitudes. Remarkably, in addition to the trivial quiescent plasma solution u — 0, there is a second 
Moo = 4-\/ (1 — cij^)/3 solution satisfying iJ(uoo) = 0. Thus, a plasma wave with H = is dynamically bi-stable: after 
the passage of the beatwave it can be either quiescent, or have the finite amplitude Uoa- It is conjectured that, by 
using a beatwave pulse with ag > acrit, the latter solution can be accessed, thereby leaving a wake of a substantial 
plasma wave with amplitude Mqo • 

This conjecture is verified by numerically integrating Eqs. (|2I3|) for two different dctunings (resonant, with uj — 1, 
and non-resonant, with uj — 0.95) and two beatwave amplitudes (sub-threshold, with oq = 0.018, and above-threshold, 
with ao — 0.023). In all cases the Gaussian pulse duration was chosen tl = 150. In physical units, for the plasma 
density of no = lO^^cm"'^ the corresponding pulse duration is = tl/uJp ~ 750 fs. Simulation results are shown in 
Fig. 121 where the solid line corresponds to the most interesting of the three cases: uj = 0.95 and ao = 0.023. The 
plasma wave amplitude of u ~ 0.75 in the wake of the laser pulse is in a good agreement with Uoo = 0.72. This 
wake owes its existence to the dynamic RB: upon interacting with the above-threshold laser beatwave, plasma wave is 
transferred from the quiescent state of m to the excited state of m = Uoo ■ The sub-threshold excitation (dot-dashed 
line) with the same detuning fails to transfer the plasma into the excited state, yielding a negligible wake that is an 
order of magnitude smaller than in the above-threshold regime. 

Linear theory also fails to describe the strong wake in this example because the detuning and the pulse duration 
are chosen such that the hnear prediction uun = ao/(l — w^) x exp [— t|(w — 1)^/4] fa is negligibly small. Resonant 
excitation (dashed line) also yields a much smaller wave. Moreover, the resonantly and the sub-threshold excited 
plasma waves would have been even smaller had the adiabatic assumption been fully satisfied. Indeed, it is numerically 
confirmed that the wake amplitudes for the resonant and the sub-threshold excitations rapidly decline for longer pulses, 
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FIG. 3: (a) Excitation of a plasma wave by a pair of identical Gaussian heatwave pulses (dot-dashed line) separated by the 
delay times Td = 920 (sohd line: wake depleted by the second pulse) and Td — 980 (dashed line: wake unperturbed by the 
second pulse). Pulse parameters: same as in Fig. |21 fflo ~ 0.023, ui — 0.95, and tl ~ 150. (b) Sequence of phase lockings and 
phase releases for Td = 920. 
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FIG. 4: Dependence of the residual plasma wave amplitude u in the wake of a pair of identical heatwave pulses on the time 
delay between the pulses td- Zero corresponds to td ~ 900/^;^. Pulse parameters are the same as in Fig.|3 



whereas the amplitude of the non-resonant above-threshold excitation is insensitive to the heatwave pulse length t^. 

To demonstrate the bi-stable nature of the relativistic plasma wake, excitation by a pair of identical heatwave pulses 
is considered. By varying the time delay Td = ujptd between the pulses, plasma wave can be either returned to the 
original quiescent state u — (Fig.|3Ia), solid line, delay time Td = 920), or brought into the excited state u = Mqo 
(Fig.EIa), dashed line, delay time Td — 980). Depending on the time delay Td, there are, essentially, only two outcomes 
for plasma wave amplitude: u « or u = Uoo- This result is remarkably nonlinear: the linear theory predicts that 
the wake behind two pulses depends on their separation in a sinusoidal way: u{t = oo) = 2u(ti)cos^ [-KTdioJ — 1)], 
where <C ti ^ Td is the instance well after the end of the first and before the beginning of the second pulse. The 
dependence of u{t = oo) on the delay time plotted in Fig. ^illustrates the effect for the identical Gaussian pulses with 
TL = 150, ao = 0.023, and w = 0.95. 

Dynamical RB described in this Letter is different from the standard equilibrium bi-stablity of a weakly-damped 
nonlinear oscillator in that the former does not require phase-locking, only the conservation of the effective 

Hamiltonian H. As Figl^^b) indicates, phase locking at 4>o — it/2 exists only during the switch-on half of the heatwave. 
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—2tl < r < 0. The plasma wave phase is released afterwards (0 < r < Td) as the pulse amplitude settles into u = Uoo- 
Depending on the delay time, the second pulse can either (i) lock the phase at = '''/2 (as shown in Fig|3|[b)), with 
the consequent decay of the plasma wave to u w 0, or (ii) fail to lock the phase, resulting in u = u^a after the pulse 
pair. Phase locking at 0o — — 7r/2 indicative of a transfer to the equilibrium Branch 3 and, therefore, equilibrium 
bi-stability, is never observed. 

So far the effect of the plasma wave on the driver has been neglected. Of course, the energy of the plasma wave 
is supplied by the heatwave. Since the plasma wave energy changes non-monotonically, different portions of the 
heatwave either lose or gain energy. In the weakly relativistic case, the plasma energy density Up ~ nQinc^v? /2. 
For concreteness, I concentrate on the above-threshold case plotted in Fig. [21(solid line). The leading portion of the 
heatwave (— oo < t < 64) contributes energy to the heatwave and is, therefore, depleted. If the heatwav e is produced 
by a laser pulse, this depletion can be described in the language of photon deceleration, or red-shifting IWilks et al.L 
In the context of the laser heatwave the red-shifting corresponds to the scattering of the photons from the higher 
frequency into the Stokes component. Assuming equal amplitude lasers, Ei = E2, the rate of the frequency shifting 
(per unit of the propagation length) can be found as —dto/dz w {tOp/ 'icojia) x d{v?)/dT. Therefore, the laser pulse is 
red (blue) shifted if du/dr > {du/dr < 0). 

If the heatwave is produced by a microbunched electron beam, the sign of du/dr can be related to the acceleration 
or deceleration gradient of the drive electron bunch Ez through 

£^WB "60 V2a(T) dr J ' 

where i?wB — mcujp/e is the non-relativistic wavebreaking electric field. Again, the sign of du/dr determines whether 
the driving bunch is accelerated or decelerated. For a microbunched electron driver consisting of femtosecond bunches 
with duration 5t <C produced by an inverse free-electron laser Srn, ^ ribo. It is estimated that in the plasma 

wave decay region of the driving bunch (64 < t < 112) the beam is decelerated at a rate of E^ ~ 30 GcV/m 
for no = lO^^cm"'^. Therefore, the marriage of the microbunched plasma wakefield accelerator and the dynamic 
relativistic bi-stability concepts yields a new advanced acceleration technique which takes advantage of the temporal 
drive beam structure to produce high energy femtosecond electron beams. 

Support for this work was provided by the US Department of Energy under Contracts No. DE-FG02-04ER54763 
and DE-FG02-03ER41228. 
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